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The present work deals with the project scheduling 
problems. Here an attempt has been made to develop a zero- 
one formulation which can accommodate a variety of practi- 
cal situations including pre-emption, concurrency and non- 
concurrency, force-finishing and variation of resource 
requirements during the erne cut ion of the activities. The 
objective functions considered are minimization of make- 
span, throughput time and penalty for late completion, v/^ 
(^Balas* and 'Resource - Time - Varying* algorithms (a va- 
riation of Balas*) have been used to solve' the formulated 
problems. The results indicate, that this formulation has 
an edge over Bowman's, and Moodie and Mandeville *s and is 
comparable with the Fischer’s and Davis* models.) 



CHAPTER 1 
INTRODUCTION 

Project Planning models based on network structure 
first appeared in the year 1958 under the acronym FEET (Pro- 
gram Evaluation and Review Technique) and CPM (Critical Path 
Method), These models have been further developed and refined 
and have been used in a variety of projects, including build- 
ing constructions, publishing ventures, research projects, 
defence works and industrial management* Both PERT and CPM 
assume the availability of unlimited resources which in many 
real problems is not a valid assumption. Hence during the 
last decade various methods have been developed for dealing 
with the case of limited resources. 

The resources of any organization can be categorised 
into five major classes. These are - manpower, machine, money 
material and (often overlooked) time. It is the interaction 
of the three main variables of time, money and manpower which 
adds intricacies to the already complex problem. We can 
accelerate an activity by using more manpower, but the excess 
availability of material may not be of any help. Under cer- 
tain circumstances, we may have to split, force-finish or 
link some activities. Apart from the above constraints we 
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may have indirect limitations that are beyond the control of 
planners hut they should, at least, receive some considerations * 
These are : weather (it affects crew efficiency, material pro- 
curement, labour turnover), space (limited working space deter- 
mines the crew size), efficiency (some crews and pieces of 
equipment are more efficient than others) and reliability 
(critical jobs should be given to more reliable resources)* 

The constrained resource problem has been categorized 
broadly into the following two types : 

1 . Maximization of some function of resource utilization 
subject to a project duration* In this case we can 
have minimization of maximum demand or the minimiza- 
tion of the sum of the squared/absolute deviation 
from the mean. 

2. Minimization of the project duration subject to a 
given resource availability. For multiproject prob- 

clems we can minimize some penalty function associated 
with the non-completion of the projects* 

The* reasons for the non-availability of the genera- 
lized analytical solution procedures for the constrained re- 
source problems are because of the difficulties encountered 
in making a mathematical formulation* Some of the difficul- 
ties are listed below: 
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Generally the activities can he started or continued 
with a reduced level of resources. This, along with 
the substitution of resources, greatly complicates 
the problem of various alternatives. 

In a network consisting of the activities of the 
various departments it may be desirable that we do 
the resource levelling for each department rather, 
than the whole network. Such problems are too com- 
plex for most of the existing methodologies . 

In most of the solution procedures we find that 
either the model grows too large for any practical 
sized problem or the solution procedure becomes very 
lengthy. 

The feasible solutions achieved from a solution pro- 
cedure may be quite off compared to the optimal so- 
lution. 

The interdependence of the activities as a result of 
sharing the same resource poses serious problems for 
the mathematical formulation. 

In case of multiproject problems different projects 
may have different priorities and this should be 
taken into consideration for project scheduling* 
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JL number of network problems may require activity 
splitting or (concurrency of activities and this has 
to be mathematically expressed. 

Keeping the above problems into consideration, an 
attempt has been made to develop a zero-one formulation. In 
this work the emphasis is on two main points : 

1. The model should be able to solve a wide range of 
network problems . 

2. The model should possess some special structure which 
may be exploited for solving any practical sized 
problem. 

A brief survey of the literature is presented in 
Chapter II. A generalized problem formulation is given in 
Chapter III. In Chapter IV, the solutions obtained using the 
proposed methodology and the discussions of the results are 
presented. 



■CHAPTER 2 
LITERATURE SURVEY 

The current approaches to solve the resource cons- 
trained network problems can be categorized broadly under the 
following headings : 

1. Heuristic methods. 

2. Combinatorial procedures; 

a. Modification of feasible schedule 

b. Construction of a schedule by means of a 
sequence of partial schedules . 

3. Mathematical programming methods. 

A brief description of these methods is given in the 
following paragraphs : V_ 

2.1. Heuristic Methods : 

( -j ) 

MAP which stands for Multiple Resource Allocation 
Procedure is a heuristic technique that considers the cons- 
traints of resource availability,, dead lines, activity rela- 
tionship etc. to determine the near optimal schedule. It is 
based on a set of eight rules which determine the crew size 
and the length of the project duration. Priorities are 
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assigned to the jots that have the same starting time accord- 
ing to the following rules : 

a) least total float, 

h) larger need for overall resource requirement, 

c) large crew size and 

d) sequence code . 

The main drawback of MAP lies in its inability to do resource 
levelling* In case of raultipro jects it attaches the same 
priority to each project. 

(o') 

J.D. Wiest’s model SPAR-1 (Scheduling Program 
for Allocation ;of Resources) handles multiple projects with 
fixed and variable crew sizes. It can handle jobs requiring 
multiple resources like machine, money, each of which are 
limited in quantity. Two cost functions used by Wiest are : 

1 . Maintaining the crews at peak loads and paying them 
irrespective of whether they are active or not. 

2. Assigning exponential increase in penalties when the 
project is not completed within the due dates. 

This model has been used in quite a number of large 
projects. It generates feasible solutions which are sometimes 
quite off the optimal one. The author has cited an example 
of this in the case of a project which involved sixteen 
# activities. It was found that the priority rule based on 
total slack resulted in a non-optimal solution. 
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(o) 

In another paper Wiest w has proposed, an extension 
of his original model. In this model (called SPAR-2) instead 
of defining a critical path, he defines critical sequence 
■which is determined hy technological orderings set of job 
times and resource constraints. Moreover, the critical sequ- 
ence depends upon a given feasible schedule. He has discussed 
several properties of the large projects. The work is exten- 
ded to the cases where resource availability is varying during 
the scheduling period. However, the author has not solved 
any practical sized problem to confirm the utility of the 
technique . 

(4) 

Burgess and Killebrew have solved the resource 

% 

allocation problem for the repetitive network. They have 
used the method of minimizing the sum of squares of the resour- 
ce requirements on various days. By moving the activities 
backward and forward (keeping in view of the technical cons- 
traints) the minimum sum is achieved. Various suggestions have 
been proposed by the author to achieve this minimum. However, 

none of these approaches guarantee convergence to optimality. 
( 7 ) 1 

Dewitte w uses the Burgess and Kilbrew’s routine with diffe- 

A 

rent measure of effectiveness; viz, the minimization of the 
sum of absolute magnitude of fluctuation from a calculated 
mean level of usage. 
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Levy, Thompson and Wiest^ 0 ' 1 have employed the Cri- 
tical Path Method of analysing the job to be done on a given 
project. Jobs are assigned at their earliest possible start 
time and the resource profile is drawn for each project. 

Then the program sets the ! trigger levels’ one unit .below 
the peak requirements in each of the projects and attempts 
to reschedule the activities. They have taken an example of 
55 activity network distributed in 4- shops. The authors claim 
that their procedure can be used for a 500 activity network. 
The effectiveness of the procedure depends upon the the amount 
of slack in the activities and the number of non-critical 
activities, N 

(*) 

Wilson has incorporated Dynamic Programming to ■ 

the Levy’s method. However, it suffers from a major drawback 

that each unit is assumed to have one resource unit only. 

(a) 

Kelley 7 has developed serial and parallel methods of solving 
the resource allocation problem. In the former, he determines 
the early start time of each activity and attempts to schedule 
it at that time. If the required resources are not available 
the activities are delayed to the earliest feasible start 
time at which resources are available. Activities can be 
:plit arbitrarily into integral time units. . Unlike serial 
methods which schedule one activity at a time, parallel 
methods schedule several activities at a time. The procedure 
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starts with the determination of the set of activities P(t), 
which are in process or could start at time t and the sub- 
set of activities Q(t) in P(t) which can he scheduled at time 
t . The author claims that serial methods are more practical 
than the parallel methods, 

RAMPS an acronym for Resource Allocation and 

Multiproject Scheduling, is a computerized method for handling 
several projects subject to some specified constraints. Three 
sets of input data are required for each activity; amount of 
resource required, time required and cost of splitting the 
activity once it has begun. Further more, starting dates of 
the projects, desired completion dates, penalty rates and 
project priorities are required. Since details of ‘RAMPS 1 
computational algorithm are not available, it is not possible 
to analyse the method. However it is claimed that it can 
solve network comprising of 700 activities with 20 different 
types of resources. 

James H. Patterson and Walter D, Huber have su- 
ggested minimum bounding/maximum bounding and search algori- 
thms to determine the duration of the project. They have 

( 19 ) 

compared their results with those of Davis 7 and Moodie 

( pQ) 

and Mandeville v J/ and they claim that their results are 


favourable . 
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2.2 Combinatorial Pro c e dure s 

a) Schedule Modifying Methods : 

In these methods we start with a schedule that satis 
fies the precedence constraints and then we modify the 
schedules so as to obtain an optimal solution which 
satisfies resource requirements. In a CPM network the 
longest path corresponds to the duration of the network 
and is computed by a simple recursive relationship. 
Similar attempts have been made for the formulation of 
an equivalent problem of finding a minimaximal path 

in a disjunctive graph. v 

( * 1-1 \ 

Bennington & Wennis v J have considered the case of 
two machine sequencing problems in which they attach 
pseudo precedence relationships for the activities 
requiring the same machine. The disjunctive graphs 
so obtained cannot be solved by CPM because the longest 
path in this case is unbounded. These problems of mach- 
ine sequencing are referred to as the two-state mini- 
maximal path problem. The project scheduling problem 
is slightly more complicated because there may be more 
than one unit of each resource available and several 
activities may overlap in the feasible solution. This 
problem is referred as a three state minimaximal path 
problem in a disjunctive graph. 
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Plorian, Trepant and McMohan^^ have suggested 
a method to solve the two state minimaximal problem of 
machine sequencing. Their algorithm for constructing 
the tree of resolution is quite effective. They first 
attempt to find a feasible solution, then attempt to mo- 
dify it and then try to prove the optimality. They claim 
that certain modifications can be incorporated to take 
into account the setup time for each operation. They 
have solved upto 10 machines - 100 jobs sequencing pro- 
blems* They claim that the solutions obtained by their 
methods are comparable to the solutions obtained by the 

Y-iq} (-Ik) 

solution procedures of Schrage f s ° and Balad.* 

Sussman^^ lias made an improvement over Plorian, 

( 12 } 

Trepant and McMohan 1 s procedure by taking the partial 

resolutions at the earlier conflicts and finding the 

( 14 } 

bounds of the resolved graphs « Balas 1 method of 
finding the resolution consists of starting with any 
feasible sequence. Each new term of the sequence is 
generated from an earlier one by complementing one dis- 
junctive arc. Here the search is drastically cut down 
because the disjunctive arcs that are to be considered 
for being complemented fall over the critical path. 

The bounds are calculated at each stage to determine 
the search direction. The author has not given any 
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example "but he claims that his method can generate a 
reasonably 'good* feasible solution. 

There is only one method, which was proposed by 

■ M < ) . 

Gorenstein v for solving the project scheduling problem 

by schedule modification. He finds that the main problem 
comes in determining the generalized coefficient of sta- 
bility GCS. for each resource j, These GCS.s are used 
not for bound determination but for identifying the re- 
source feasibility of the resolution. Further more, if 
GCS. is less than resource availability it implies 

d • 

resource feasibility, but the converse is not necessa- 
rily true. ■ Hence it may be possible that the optimal 
solution is not achieved in using the GCS . criterion, 

d 

b) Schedule Construction Methods : 

As the name implies , the general approach in these 
algorithms is to 'build up* a feasible solution. In 
almost all of the following methods an ordering or 
permutation of the jobs is tried to build up the sche- 
dules, Schrage ~ >J has evolved a branch and bound 
method for implicitly enumerating all schedules and 
determining the optimum. He has considered the follow- 
ing constraints: (1) Precedence (2) Resource availa- 
bility (3) Won-preemption. 
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The procedure suggested is essentially a tree gene- 
rating procedure. For each partial schedule or branch in 
the tree one can calculate lower bounds for all complete 
schedules that can be generated on that time. For obtain- 
ing the bounds on the tree he suggests (1) Precedence 
based bounds, or, (2) Resource based bounds. The objec- 
tive function is min. ( max. finish time) . The gene- 
ralizations are applicable to the following cases 
(a) multire sources (b) varying resource availability 
(c) ’OR 1 activities (d) unavailability of resources for 
changeover between activities (e) wide range of objec- 
tive functions . 

(17) 

Schrage has suggested improvements for taking 
preemptive cases. For the objective of min. (max. finish 
time) the bounds are either precedence or resource based, 
Ee indicates that the preemptive version (a) requires 
slightly more computation (b) does not result in a subs- 
tantially lower min-max completion time. Mason and 
Hoodie l s v J method minimizes the combined cost of flu- 
ctuations in resource demand and delay of project dura- 
tion. Resource requirement and duration of each activity 
is represented by a vector. At each time the amount ir >f 
work remaining and the total cost incurred are deter- 
mined. These costs are used for further branching of 
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the tree. The resource availability constraints can he 
used at a node to eliminate some subsets of the tree. 
According to the authors, the results indicate that the 
optimal schedule is relatively insensitive to changes in 
the ratio of project delay cost and the resource fluc- 
tuation cost* 

Davis and Heidorn's^ ^ algorithm permits the com- 
putation of optimum project duration under the conditions 
of multiple resource requirement per job. The sample 
network is converted into unit duration activity network 
by imposing the constraints *must immediately precede 1 
along with * normal precedence*. The number of stages 
are equal to the number of time periods in the longest 
path. Feasible subsets are generated and from them r A 
network* of the sample problem is constructed. Direct 
arcs are drawn between nodes if precedence and resource 
constraints are met. The trouble arises due to the 
fact that if activity durations are large the number of 
feasible solutions shall be extremely large and hence 
the identification of ’A-network* will be problematic. 

The authors have suggested the use of dynamic progra- 
mming approach for making * A-network* from the feasible 
subsets. Branching is simplified by using, (1) pre- 
cedence based or (2) resource based criteria. This 



algorithm can also account for the -variations in jot 
continuity and varying resource requirements. 

( 20 ) & 

Johnson has formulated branch • bound algorithm 

A 

for the construction of the schedule. Johnson 1 s algorithm 
uses the following logic for the construction of the 
schedule : 

1 . Given a partial schedule at time t^ consider 
the set of all /unschedule activities whose 
immediate predecessor have been scheduled. 

2. From (1) form subsets of activities so that 
resource constraints are considered. 

3. By including these subsets form new subsets. 

He considers two bounds for each partial schedule, 

(a) the critical path gives one bound and, (b) resource 
requirement of the other. Based on computational 
experience, Johnson claims that his is the most com- 
prehensive algorithm for scheduling single resource 
problems • 

2.3 Mathematical Programming Methods i 

( 22 ) 

H.H. Greenberg has developed mixed integer for- 
mulation for the general n job - m machine sequencing prob- 
lem. He has used S . and d. to denote the start time and 

xp xp 
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duration of jot i on machine p. Further he defines t- 

r 1, if job i starts on machine p before job j. 

! 1 

0, otherwise. 

With these he has specified the technical constraints and 
machine availability constraints. By dropping opt the inte- 
grality restrictions the problem is solved by branch and bound. 
It is worth noting that the linear program solved at each node 
is nothing more than a CPM problem for which we have better 
methods in disjunctive graph technique. 

Bowman (24-) gj_ ven the zero-one formulation of 

schedule sequencing problem. He has considered the example 
of 3 products which have to be processed on 4 machine in a 
given technological ordering. He uses = 1 to denote 

that a product x is being processed on machine A at time i 
and X A , . = 0 otherwise. A large penalty is attached to the 
jobs being processed at the end and this helps in finding the 
shorter make span of the project. The number" of variables 
used in his method are extremely large (are equal to 
(products) x (machines) x (time periods)). The details of 
this method are given in appendix A. 

Marine has suggested an improvement over this 
Integer valued unknowns 1 s are used to indicate 



method. 
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the day on which task j is to he started. He has converted 
the noninterference constraints into the equivalent constraints 
by introducing the zero-one variable lb, for jobs j and k. 

jK 

These modifications are used for expressing the sequencing 

constraints, due date requirement* and exact delays between 

two consecutive jobs . The number of variables required in 

( 24 -) 

Marine’s method are considerably less than in Bowman’s 
method. Manne has suggested that all attempts should be made 
to reduce because these are required only in noninter- 

ference constraints and in any numerical problem these classes 
may occur only in a few cases. In Wagner's method, the 
number of constraints to be satisfied are staggering but most 
of them are inoperative.. He claims (without any proof) that 
his model can solve a wide class of machine sequencing prob- 


3ms . 


( 24 ) 

Hoodie and Hands ville ~ >J have solved the resource- 

( 27 ) 

balancing problem by using the Bowman* s v ' ' method for solv- 
ing the assembly line balancing problem. The difference 
between the two class of problems is that an activity in 
assembly line balancing problem occurs at one and only one 
station while, in resource balancing problem, an activity may 
spread over several days. An Integer Linear Programming 
formulation is used to minimize the cost of applying the 
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resources* Sven the multiproject- resource-balancing' problems 
have been solved by these methods. The details of the for- 

('po> 

mulation of Hoodie and Mandeville are given in Appendix E. 

/pO\ 

Fischer^ J has used Lagrange multipliers to obtain 
a strong lover bound on the cost of an optimal solution for the 
resource constrained network scheduling problems. A branch 
and bound procedure is used to obtain the lower bounds* The 
concept of' active schedule has been used to eliminate nodes 
in the tree and to update multipliers at each branching. The 
author has suggested that the efficiency can be improved by 
the elimination of inactive scheduling. He assumes that one 
type of resource is available at each time period. The 
details of his method are given in Appendix C. 



CHAPTER 3 


PROBLEM FORMULATION AMD SOLUTION METHODOLOGY 

The generalized problem formulation of a project 
scheduling problem is presented in this Chapter. 

3.1 General Statement of the Problem ; 

! Eor a project scheduling problem, minimize either 
throughput time or make-span or penalty for late completion 
subject to the various constraints. The various constraints 
considered account for the concurrency and non-concurrency of 
activities, availability of resources, technological ordering 
preemption, force -finishing of activities, etc.’ 

3.2 Nomenclature : 

The following nomenclature is used to define the 

problem ; 

i : project number ; i = 1,2,..., I. 
j ; activity number; j = 1,2,..., ; 

KB = number of activities in project i. 

t : time period. 

T : period for completion of all the projects. 

D . . : duration time of activity j of project i. 
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EA, 


iD 


M.. 


EP ± 


LP. 


earliest finish time of activity j of 
project i. 

latest finish time of activity j of 
project i, 

earliest finish time of project i. 
latest finish time of project i. 


^ijkt 1 * : amoun ^ resource of type k required ty 
activity j of project i at time t l . Here 
t l is the time since the activity has begun, 


resource number, k = 1, 


• > 


K; where K is 


"it 


T. 

ic 




the number of different types of resources 
required. 

amount of resource of type k available at . 
time t. 

cycle time after which the project i is 
repeated. 

available which takes the value 0 or 1; 

X^ = 1 indicates that the jth activity 
of project i is completed in time t and 

is zero otherwise. X^^ is a variable 
only from EA . . to LA. . because it can be 

J-J -LJ 

assigned a value of zero for t < EA^.. an( 3- 
t > LA ± ^ . 

"If the resource requirement of an activity does not vary 
during its duration, we shall use only 
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Ht i a variable -which takes the value 0 or 1 ; 

1^+ = 1 indicates that the project i is 
completed by time (t - 1 ) and zero otherwise', 
^it is a variable from (EP^+1 ) to LP^ only 
because it can be assigned the value zero 
for t < EP_^ and one for t > LP^. 

Z t : a -variable which takes the value 0 or j; 

Z^=1 indicates that all projects have been 
completed by time (t - 1) and zero otherwise. 
It is a variable for t > max EP^ to t ^ T . 

Pit : penalty incurred due to non-completion of 
project i in time t. 

r ■ 

3 ,3 Assumptions ; 

The following assumptions have been made for the 
mathematical formulation of the project scheduling problems 

1 • The time duration of each activity is 
deterministic, 

2, No crashing of activities is allowed, 

3* Resource requirement of each activity is 
known, 

4-, Resource availability throughout the project 
duration is known. 
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5* Duration time of each activity is independent 
of the other activity, 

6, No alternate sequencing of an activitjr is 
permitted, 

7* All activities in a project are equally impor~ 
tant , 

8, Resource substitution of an activity, is 
not allowed. 


3.4- Constraints : 


The various constraints can be expressed mathe- 
matically as follows : 

1 • Activity Completion Constraints 

Each activity of the project has to be completed. 
Mathematically this can be represented as. 


LA. . 

t=EA. . 
id 


Then, by definition 


■^IjjLAij 


XA. r 1 

T 4 , 


t=EA, 


X ijt 


n 


(D 


( 2 ) 


Since can be either 0 or 1 , we can xvrite Eq.(1) 

in an equivalent form as given below: 


23 


'LL. .-1 
3-3 

^ X, . 4 . < 1 

t=EA. . ^ 

3-3 


( 3 ) 


The formulation as given in : Eq, (3) has the advantage 
that the number of variables in an activity is reduced 
by 1 . 


2* Project Completion Constraints : 

for the completion of project i by time <t - 1 ) 
we must complete all its activities in time (t- 1 ) . 


Mathematically 


N. t ~ 1 

N. 1.4.4: X- - 

' L li: 3=1 P =EA iij 1 ‘ )p 


i — 1 j • • » j I 3 


( t > EPj 


t^ T 


Gf) 


If all the projects have to-be completed by time 



we. re quire f 



1 

JLs 

x . 4 i, 

.3. i=1 

TT 

% 

Jifcj 

3rg 

t=1 

3=1 

p=EA. 


3-3P 


t > max 
t 4' T 


( 5 ) 
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3» Force-Finish Constraints ; 


An activity ij is said to be force -finished if 
it has to be force finished by time t^ (EA^ < t^.,<. 
LA. .). To satisfy this constraint we modify Eq. ( 1 ) 

-i_J 

as follows : 


*f 

I 1 . 

t=EA. 


x ±dt 


10 


( 6 ) 


4-.. Sequencing Constraints ; 


Let activities 3 ^ and j^ ^e the ^ wo activities of 
the project i. If j.j precedes j 2 , this can be repre- 
sented by the following inequality: • 


1 * x i 3 t + D i 3 , 

i 


“« 2 




t=EA. 


tX ±3 2 t 


10 - 


(7) 


The precedence relationships between the activities 
of the different projects can be expressed as, 


LA. . 
10 

T! 


t X 


t=EA. . 

-L J 


i-jO^’t 


+ D. 


L &4 4 H- 


1 2 3 2 ' 


t X. 


t— Eil . • 

3-2 3 2 ^ 


i 2 5 2 t 


( 8 ) 


where i^ and i^ refer to different projects. 
Furthermore, if we want to have some lag (sayV) 


between the two consecutive activities we can express 
it by imposing the equality constraint. Mathematically 
this is expressed as, 


LA. . 

2 1 


LA 


i 2 3 2 


t X -‘ ii + D T.i. + v- 


i-i 3 




X 2 3 2 


t=EA. 


(9) 


3 - 2 3 2 ~ 2' 1 2 


5*. Concurrency and Non-Concurrency Constraints : 


Two activities are said to be concurrent if both 
of them have to be carriedout simaltaneously. If . eithe 
of them has to be completed before the other starts 
we call them as non- concurrent activities. To ensure 
that an activity ij is being carriedout at time t we 
must satisfy the following relationship. 


t+D. .-1 

J- J 

n 

p=t 


x. . 


= 1 


• * ♦ • 


(10) 


Hence if we want that only one of the activities i-jO-j 
and i- 2 J 2 be carriedout at time t # we" should satisfy 
the following inequality. 


2 

p=t 


“1 

3-jP 


t+D. . -1 
^ 2 J 2 

* ■ 

p=t 


XjL 2 3 *2 p 





on 



j 


r-v /* 


For non-concurrent activities Eq. (11) should he satis- 
fied for the following range of t. 


Min. (LA.. . LA. . ) > t Max. ((EA- • - B . • + 1), 

■J-jj-p ' l 2?2 i j 1 - U [ J -] 


(EA. * - B . . + 1)) 


^ 2^2 1 2‘*2 


For concurrency we need to add the resource requirement 
of both the activities. Mathematically it is expressed 
as given in Eq, (12). 

6. Preemption Constraints : 

Preemption of an activity implies that an activity 

can he interrupted, once it has started. Let us say 

that an activity ij of duration B . . is interrupted 

after duration B. . and then restarted and completed 

without any further interruption. Thus the activity 

i j is divided into two subactivities i and i^ °E 

durations B. . and B . . respectively, such that 

D. . = D. . + B.. . On these two activities wecan 

13 13-j 13 2 

easily impose the sequencing constraints as in (8). 
Similarly we can take into account more interruptions 
of an activity. 
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7. Resource Availability Constraints : 

If the resource requirements of an activity do 
not vary during the execution period of the activity, 
the resource requirements at any time t wopld he given 
t>y, 

I . H. t+D^-l 

B % hi,* “it * = 1 ^ . 


k 1 , • * * , K 

If the resource requirement of an activity changes 
during its execution period, we can satisfy the 
resource requirement at time t as follows : 


( 12 ) 


I 

% 

t+D ij 

21 

% 

21 

i=1 

1=1 

P^t 


X^R/L, t = 1 j • « • jT ; 


R ijkt' X ijp- M kt ‘ } 


(13) 


k = 1 , . . . ,K 

For the repetitive project the resource requirement 
has to he satisfied only for 1 time units. 

3-7 Objective Functions ; 

The various objective functions considered in 
this study are concerned with the minimization of make span, . 
throughput time and penalty. Ma thematic ally they can be ex- 
pressed as follows : 
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1 , Minimum Make span: 


The completion of all the projects in the minimum 

time can he expressed by the following objective 

j ‘ 

function : 


i 

Maximize Z = IE) Z±. (14-) 

maxJ2P^+1 


2* Minimum Through-put- time : 


The earliest completion of a project i can he 
achieved hy the following objective function: 

' LP i 

Maximize Z = 2D ^it (15) 

t=EPi+1 


The throughput-time of all projects can he minimized 
hy the following relationship : 


I LP ± 

L&ximize Z - Y it 

i=1 t=EP ± +1 


(16) 


3* Minimum Penalty : 

/. 0 

Let P^j. he the penalty for the non-completion of 
project i in time t, then we can express the minimum 
penalty as follows : 

1 LV 1 

Maximize Z '=** ^ P.**. Xj-t. (17) 

i=1 t=EP.+1 it 
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3*6 Solution Methodology : 

Using the aforesaid generalized formulation a 
number of project scheduling problems have been formulated. 
The descriptions of these problems are presented in Chapt- 
er IV. The following algorithms have been used to solve 
the formulated problems, 

1* Balas* algorithm for zero-one programming has been 
used to achieve the optimal solution. Since it 
takes many iterations to find the optimal solution, 
another algorithm is proposed. 

2. Resource - Time - Varying algorithm, a variation of 
Balas 1 algorithm has been developed by which we gene 
rate the feasible solutions in such a way that they 
are optimal. The various steps required are given 
/" below: 

a) Schedule the projects based on the early start 
times and find the maximum resource requirement. 

b) Decrease the resource availability by one and 
find a feasible solution. If a solution is ob- 
tained, decrease the resource availability furth- 
er by one and continue till WSL. . don’t get a 
feasible solution. 
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e) Increase the completion time of the project in 
step of 1, till a feasible solution is obtained, 

d) Repeat steps (b) to steps (d) until the resource 
availability is equal to the maximum resource re- 
quirement of an activity or when the resource 
utilization has reached a reasonable value . 

The advantages of using the second algorithm are 
(1) we determine a schedule for each value of resource 
availability f (2) at each step we are finding a feasible 
solution (which is optimal) and this is computationally 
faster than finding the optimal from the feasible solutions, 
and (3) in step (b) we are doing the resource levelling. 



CHAPTER If 


ILLUSTRATIVE PROBLEMS 

In this chapter a few problems have been solved 
to illustrate the proposed formulation. The networks have 
been taken from various books and journals. For each prob~ 
lem, the network, the scheduled graph and the resource 
profiles have been developed. 

Network 1 : 

An 11 activity, 1 resource, activity on node net- 
work is given in Fig. 4*01. The problem has been solved by 
the Resource-Time-Varying algorithm. The results are given 
below : 

The scheduled graph and resource profile based, on- 
early start time are given in Fig. *+.02. The maximum re- 
source requirement is 9, The resource availability is 
decreased to the level of 8 and a feasible solution is sea- 
rchod for T = 18. No solution exists for T = 18 and so we 
increase it to 19 and obtain a solution as shown in Fig 0 
*+. 03. For this schedule the maximum resource requirement 
comes out to be 7 only. Further, on decreasing the 




<e source availability to the level of 6 we get a feasible 
solution at T = 21. This is shown graphically in Fig. 4.04-, 
The total number of iterations required for the solutions 
are 1270. The resource utilization for the three solutions 
are 111/162, 111/133 and 111/126 respectively. | 
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5 

, ? 7 > ’—7 » 

i D 

IS 

1 8 Time 



Resource Profile 

Fig* .02 













36 



Network 2: 

A 9 activity, 2 resources, event-on-node network 
is given in Fig, 4.0^. This problem illustrates the 
minimization of make-span and throughput times for the 
case of non-re petitive and repetitive projects respectively 
For the resource availability of resource type 1 and 2 as 
3 and 1 respectively, the scheduled graph and the resource 
profiles are given in Fig. 4-. 06, The make-span obtained 
is 12 units of time and the nufnber of iterations required 
are 1270. When the project is repeated every *> units of 




time, we get, for the resource availability of type 1 
and 2 as 6 and 3 respectively, the. scheduled graph and 
the resource profiles as shown in Fig* 1.07 . The through- 
put time is 15 units and the number of iterations required 
as?e 373a The resource utilization of type 1 and 2 are 1/5 
and 8/15* On decreasing the resource availability of type 
2 to 2 (because its utilization is lower) we don't get any 
feasible solution for the maximum allowable throughput time 
of 15 units o Next, taking the resource availability of type 
1 and 2 as 5 and 3 respective!}^ we get the scheduled graph 
and resource profiles as shown in Fig, 1.08, The through- 
put time is 15 units and the number of iterations required 
are 1688 e The resource utilizations for resource type 1 
and 2 are 21/25 and 8/15 respectively. 

Network 3 : 

A 2 Jobs ,3 machines flow-shop problem is given in 
Fig. 1*9 • The total completion time of Job 1 and Job 2 
are ll and 12 units respectively* The penalties incurred 
for the non-completion of Jobs in time t are as given 


below ; 
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Resource Profile 2 
Pig. V.06 

Using Balas* algorithm for minimizing the total penalty 
we get the scheduled graph as shown in Pig* The 

number of iterations required are 31 J> 9 * 




^Network starting 

r 


f Network starting 
V 5 time units before 


7 Network starting 
1 10 time units before 





bo 


* i rty<r 



Scheduled graph 


Network startint 
at t = 0 


Network starting 
J? units before 

Network starting 
10 units before 


l 



Resource Profile 1 



•h lime 
Resource Profile 2 


Fig. 4.08 



Network 4 - i 

A 3 jobs, 3 machines job-shop problem is given in 
Fig c 4- *11 * The sequences of the jobs on the machines are 
pre specified as follows : 

Machine Sequence for job 1 i (1, 3, 2) 

Machine sequence for job 2 : (3 , 2, 1) , 

Machine sequence for job 3 s (2, 1,3) 

Minimizing the make -span of the project we get 

the scheduled graph as shown in Fig. 4-. 12. The make-span 
is 13 units of time and the number of iterations required 
are 5*6. Both networks 3 and 4- assume availability of 
1 unit of resource (machine) only. 

Network 5 : 

A 22 activity, 1 resource event-on-node network 
is given in Fig. 4-. 13* For the project duration of 30 units 
the attempt has been to achieve the maximum resource utili- 
zation. Decreasing the resource requirement from 9 units 
(corresponding to early start schedule) to 6 unit we . 
obtain the scheduled graph and the resource profile as 
given in Fig. 4-.14-. The number of iterations required 



4*2 


Machine A Machine B Machine © 



5 3 




Joh 1 

Job 2 

Job 3 
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Network 6 i 

An example of inter related project requiring two 
resources is given in Fig. 4.15* Projects 1 and 2 consist 
of 7 and 9 activities* Activity 2 - 4 of project 1 has 
to precede activity 10 - 12 of project 2. For the resource 
profiles as shown inFig, 4.16* The number of iterations 
required are 5016. 

Network 7 : 

A 16 activity, 1 resource network is given in 
Fig. 4.17® In this example activities 1-2, 3-8, 5-7 
6-9 and 8-11 have to be carried out in the first 
department , and the remaining in the other* The resource 
availability for the first department is 2 units while 
for the second department it is 3 units. The scheduled 
graph and the resource profile for the minimum make-span 
are given in Fig. 4,18. The number of iterationsre quire d 
are 81 • 
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CHAPTER 5 

INFERENCES AND SCOPE FOR FURTHER WORK 

The following inferences have been drawn from the 
study of the generalized formulation and the solutions of 
the project scheduling problems. 

1. The Balas* algorithm, used in this work, turns out 
to be quite efficient for finding the feasible so- 
lutions but it takes many iterations to obtain op- 
timal from the generated feasible solution. 

2, The proposed formulation requires (HI HL 

•V"! -VjJ 

(I4L. - EA. .) + y_ (LP. - EP. )) variables for the 
project sche dialing problems for minimum penalty 
and throughput times. The nunber of variables, 
when the objective is to minimize make-span are 

- EA. .) 4- (T - max.- EP.)). However, 

i J . -**• ' , 

and, Hoodie and Mandeville l s formula- 
tions require (I , 211 N i • ■v r abiables. Thus, 
the proposed algorithm requires much less number 
of variables than required by eithpr of the two 
algorithms. Therefore this algorithm can be used 


c n z_ cia.. 



to handle larger problems efficiently « However s 
as the number of activities increase beyond a cer~ 
tain limit the model becomes unmanageable,. Furth- 
ermore it should be no.ted that when the duration 
of the project is increased by one unit of time,, 
the number of variables increase by ( lb + lj 
for the objective criteria of minimum penalty and 
throughput time. In case of minimum make-span tnoy 

increase by (T1 N., + 1)* 

i=1 1 


The proposed formulation has an edge over 
Fischer* formulation. This is because in 
the latter we have the assumption of availability 
of one unit of resource throughout the project 
duration. When this assumption is relaxed the 
efficiency of his method goes down while it is not 


so for the proposed algorithm. Furthermore, 
Fischer’s model is useful when there are many pro- 
jects with relatively less number of variables in 
each of them while’ the proposed formulation is 
preferable when the number of activities m a 


project are large and there are fewer projects. 

This is due to the fact that in Fischer’s method 
the number of feasible schedules iip|^e^r^__. 


Abe. No. A 42654] 
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rapidly with the increase in number of activities 
in a project. 


If. A comparison of the proposed algorithms with Davis 
and Heidorn T s 

algorithm indicates that the 
efficiency of the former does not go down for the 
large activity networks, as is the case with the 
latter. 


5# The Resource - Time - Varying algorithm suggested 
in this dissertation has its relevance for the 
specialized 2 *• dimensional cutting stock problems 
(e.g., wood; where the grain structure forces the 
sawing to be done in a special direction only). 

If the total length of the plate, total width of the 
plate, the length of the plates to be cut, the 
•width of the plates to be cut and the order in 
which they have to be cut are to be designated by 
project span, resource availability, resource 
requirement, activity duration and sequencing 
constraints respectively, then the proposed algo- 
rithms give the minimum area of the plate. 

As illustrated in Chapter IV, the generalized 
formulation proposed in this dissertation can 


6. 
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accommodate a variety of project scheduling 
problems , 

cone for Further Work : 

The proposed algorithm can he further improved by 
ncorporating the following features : 

1 . Activity completion and sequencing constraints of 
the proposed formulation possess a special struc- 
ture. Furthermore, for all the projects (except 
for inter-related projects), the activity comple- 
tion and sequencing constraints are independent 
while the project completion and resource require- 
ment constraints are dependents. These facts can ■ 

be used to develop efficient algorithms, 

■ ' : ■■J 

' '■ ' \ ■■■ . \ ■ ■' ■ - " •- ; : - ' , ■ - v.V v ,, ; 

' 1 

2* The flow-shop and the job-shop problems possess 
a special structure that, except for the first 
and last activities, all other activities have 
exactly one predecessor and one successor. Further- 
more the resource availability, is 1 unit through- 
out the project duration. These properties should 
be utilized in accomplishing faster computations. j 

3, The proposed formulation yield a very sparse 

matrix. Hence the properties of sparse matrices 



should be used to achieve computational efficie- 
ncy. 

Certain modifications should be incorporated to 
take into account the resource substitution and 
the crashing of the activities. 
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APPENDIX A 



Consider the case of three products X, Y, 2 which, 
are to be processed on fopr machines A, B, C, D, The se- 
quences for products X,Y,Z are A B '-*• C -» D, C — * A~> 
D -+ B and D — r A respectively. The project is to be 
completed by time T, Let X^ , t = 1 denote that product X 
is being processed on machine A at time t and let X^.^= 0 
otherwise. The various constraints are formulated as 
given below; 



If d^ denotes the duration of the product X on 
machine A, then the job completion constraint is ex- 
pressed by the relationship; 

X, cLr. for all products and machines, (1) 
1=1 ,;L ju *’ 

2, Non-interference Constraints : 

This implies that two or more products cannot be 
handled simaltaneously on a machine. Expressing this 
mathematically we get, 
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^ *■£ ^ for t = 1 j*«*jT and all machines , 

3* Precedence Constraints : 

To ensure that product X is processed first on 
machine A and then on B we must satisfy the following 
inequality. 

3-1 

d XA 2* X*. for j = 1,...,T and all sequences? 

1=1 

For products Y and Z the precedence relationships can 
he expressed similarly. 

A. Non-interruption Constraints : 

The jobs, once started, should not be interrupted. 
Mathematically this is expressed as, 

d XA X Ati * d Xil X A:i+1 + X A:j 4 ^ 

for j = 1, *#.,T and all jobs. 

Objective Function : 

For minimizing the make span of the project the 
following objective function is used. 

Min. O.F. a 1 (Xp j23 + y b:23 + Z A23^ + ^ ^ X D :2A + X B:2 ) + + 

Z A . 2lf ) + 16 i%. 2 f Z a . 25 ) +...+ KjCXd;! + 

Y B:T + Z A :T^ * 
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APPENDIX B 

MODDIE AND MNDEVILLE * S FORMULATION 
FOR PROJECT RESQURGOSiJflCING 

Moodie and Mandeville have used the following 
notations for resource balancing problem. 

Let denote resource level of activity i on 'day 
A; resource level of activity i on any day; ,6^ duration 
of activity 1 in days and ^ maximum resource available. 

The various constraints are expressed as follows : 

1. Activity Completion Constraints ; 

Each activity must be completed in the specified 
duration. To ensure this we require the following 
rela tionships : 

A^ + + C^ + . . . + G.j ~fi ± 'T i 
Ag + Bg + C2 + . • . + Gg = ^2 ^ 2 

2. Resource Constraints : 

Each activity should be performed at the specified 
level of resource requirement. Ite. thematically it is 
expressed as, 
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A 1 + ^ ~ 1 9 ^T 1 B 1 + I B 1 = 1 » t 2 + J A 2 = 1 * 

where I’s are either 0 or 1 . 

3. Precedence Constraints : 

The precedence relationship can he represented 

as follows : 

1 B ^ - 1 * JL» 

T 1 2 ^ 1 

c 2 < <1^ A i + 'f7F 1 B i 

A* Re sou? ce Balancing Constraints : 

To balance the resource level over the total 
project?; we m$st have, 


A 1 


A 2 


+ Aq 


0 



B 1 

+ 

b 2 

^ ^ # 

+ Bq 

w** 

Aj , 

+ A 2 

+ **♦ + 

c 1 


°2 

» 

♦ 

+ Cg 


B 1 

+ b 2 

4* # * • •+ Bi 

0 

: 

% 

+ Ng 1 

+ ••• 

+ 1 

'8 

' ’ V ; . 



5 « Hon~i nte rruotlon Constraints : 

An activity once started should he completed 
without interruption* To ensure this we have, 

( h 1 ~ 1 ) ^ ^ B 1 j ( fi 1 - 1 ) B 1 $. G 1 + A- 
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Objective Function: 

The objective function used by him is concerned 
with the minimization of cost of applying the 


resources • 



APPENDIX C 


FISCHER.* S FORMULATION FOR OPTIMAL SCHEDULING 
~ USING LAGRANGE MULTIPLIERS " : 


Fischer has developed a zero-one formulation for 
the project scheduling problems using Lagrange Multipliers, 
Various notations used b& him are, 


i 

t 

k 

n i 

p ij 

%j 

f ij 

C i1 


: number of projects; i = 1,„*.,I 
: time period ; t = 1 , T 
: types of resources; k = 1, ...,K 
: number of activities on project i 
: duration of activity ij 
; start time of activity ij 
: finish time of activity ij 
: cost associated with qifch job-schedule of 
project 1 

: the amount of resource type k required at time 

. ■ ■, :■ /.■ ■ v 

t for qth schedule of project i. 


The various constraints are expressed as ..follows : 

1. Non-preemution Constraints : An activity,, once 
started, should not be interrupt^d.This can be 


expressed as follows : 
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'Ij 


*13 + 


(1) 


2. Precedence Constraints ; 

If an activity ij follows the activity ii, it 
mast satisfy the following inequality. 


t ij ^ 


'll 


^il + p il 


( 2 ) 


3* Project Completion Constraints : 

If f^ = max > then we have 


3=1,. 


j l • * j 21-^ 


f ± '$ I 


(3) 


Let a job schedule for the project 1 be any schedule 
that satisfies the constraints (1), (2) and (3)« 

Let X i = 1 Indicate that we use the qth -schedule for 

xq 

project i. 

* 0, otherwise. 

Then we have. 


Minimize £ C iq 


s.t, 




Iq 

•mA' 

► 


i 


x. = 1 i - 

^ «A»vjL . 




(A) 


^ * = V*”» T > 15 - V-v 1 


X^ = 0 OP 1 q - 1,...% > i 55 
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The above problem can be converted into the following 
Lagrangian problems. 


W(U) = fin. ( C. X- + a U, ,( TW. ■ . P )) 

X iq X i^ q kt kt \t ;; 


s.t. ^X ±q =1 


1 


X iq = ° ° r 1 


a 1 , « » . , 1 


q = 

i = 1 , * . *,I 




(B) 


Now, let 7T iq (0) = c lq + t cC. iqkt H kt and 


IT. (U) = Min. TT . (U) 

q j-q 


We can show that the dual of (A) with integrality 
relaxed is, t 


Maximize W = 5j]fT - 2 E.. U~. 

i kt K Jv 

s.t. TT i^r C iq + ^ iqkt U kt q=1 ’"’ Q l 

U kt> 0 


i— 1 j • • • j I 


k = 

t ~ 1 , • . ,T 


(C) 


J 


Where 7^, i = 1,...I are unrestricted in sign. 

: Lagrange Multipliers are positive. 

Even 'though (C) has numerous constraints,' most of 
them are non-binding at a particular extreme point 
solution* This fact is. used in developing a branch 
and bound algorithm in which we solve (C) by simplex 
method dealing with only a subset of binding Constraints 
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